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Description 

METHOD OF SYSTEMS ANALYSIS FOR DETERMINING MEASURED VARIABLES 

5 The invention relates to a method and a device for designing a 
technical system, and a corresponding computer program product. 

In the. description of a technical system, for example a power 
station, reference is made to various parameters such as,' for 

10 example, pressures, mass flows, etc. The parameters obey 

certain physical laws, such as for example those of mass or 
energy balances, which can be expressed by a system of 
equations. The solutions of this system of equations are the 
state variables of the technical system. From these state 

15 variables it is possible in turn to calculate diagnostic 

variables relevant for the operation of the technical system, 
such as for example the efficiency of a power station. The 
specific state of a technical system can, furthermore, be 
sensed by measurements. The measured variables of the 

20 measurements could directly reflect the value of a state 

variable; but it is also possible for measured variables which 
are derived from the state variables to be measured. For 
example, it is possible to measure the temperature of a 
technical system whereas the actual system state variables are 

25 the enthalpy and the pressure. In order to determine the state 
variables from the measured variables, one generally carries 
put a measurement and looks for the values of^the state 
variables which solve the system of equations, -and for which 
the derived measured variables lie closest to the measured 

30 values determined by the measurement (see for example VDI 
Guideline 2048) . * 



The problem can arise, due to the number of equations in the 
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system of equations being too small, or too small a number of 
measurement points, that individual state variables or 
individual diagnostic variables remain indeterminate . In 
addition, the state variables or diagnostic variables, as 
5 applicable, may. be subject to great uncertainty, because of 
measurement errors. It is therefore necessary to decide which 
measurements would permit. the precision of certain state 
variables to be improved, or would permit any determination at 
all of certain state variables. For this purpose, it is usual 
10 to fall back on the advice of experienced engineers, and the 
suggestions of these engineers can be checked by simulation 
programs. However, this requires time-consuming analyses. 

The object of the invention is therefore to specify a method 
15 for designing a- technical system by which a systematic 

determination is made of how the measurements of individual 
measured variables influence the parameters of the technical 
system. 

20 This object is achieved in accordance with the features of the 
independent claims. Developments of the invention are defined 
in the subclaims. 

The method in accordance with the invention is for use in 
25 designing a technical system which is characterized by 

parameters, comprising state variables and diagnostic variables 
which are dependent on the state variables. Here, the term 
designing means, in particular, the analysis and/or changing of 
the technical system, in particular the analysis and changing 
30 of the measurements made in the technical system. In this 

situation, the technical system is specified by a system of 
equations, whereby the state variables are the solution of the 
system of equations. Incorporated into the design of the 
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technical system is a measurement 'park' which includes first 
measured variables, whereby these first measured variables are 
measured in the technical system with a prescribed precision. 
In addition, second measured variables, which are dependent on 
5 the state variables, can be measured in the technical system 
with a prescribed accuracy. 

In the method in accordance with the invention, first 
sensitivity variables are determined for the first measured 

10 variables and/or second sensitivity variables are determined 
for the second measured variables whereby, for the purpose of 
defining the first sensitivity variables, a determination is 
made of the extent to which a change in the accuracy of the 
measurement of the first measured variables influences at least 

15 one selected parameter, and for the purpose of defining the 
second sensitivity variables, a determination is made of the 
extent to which the measurement of the second measured 
variables influences at least one selected parameter. The 
measurement park is then amended, depending on the first and/or 

20 the second sensitivity variables, in such a way that the 

accuracy of one or more of the first measured variables is 
changed, and/or one or more of the first measured variables is 
taken out of the measurement park and/or one or more of the 
second measured variables is added into the measurement park. 

25 This amended measurement park is used for designing the 
technical system. 

In a preferred form of embodiment, the accuracy of a first 
measured variable is preferably increased, if the first 
30 sensitivity variable for this measured variable lies within a 
predefined value range, and/or a first measured variable is 
taken out of the measurement park if the first sensitivity 
variable for this measured variable lies within a predefined 
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value range and/or a second measured variable is added into the 
measurement park if the second sensitivity variable for this 
measured variable lies within a predefined value range. It is 
thus possible in a simple manner, by the choice of different 
5 value ranges, to modify the design procedure appropriately for 
various user-specific requirements. 

In a preferred form of embodiment, the technical system is 
specified by a system of equations H{x) = (Hl(x) ... , Hn(x) = 
10 0, where x = (xl, ...xn) is a vector which includes as 

components the state variables xi . It is noted at this point 
that all the indices i, j, k or 1 used below represent cardinal 
niimbers . 

15 For the purpose of carrying out the method in accordance with 
the invention in a preferred form of embodiment the following 
matrices are, in particular, calculated: 

a matrix N which spans the null space of the Jacobian matrix 
20 H, 

a matrix W, such that Iff-W is the inverse of the covariance 
matrix of the first measured variables yi=bi (x) , where the 
covariance matrix has as its entries the covariances 
aij^=E((yi - E(yi))(yj - E(yj))), where E(y) is the expected 
25 value of y; 

a matrix M which is the pseudoinverse matrix of A==W'Db'N, 
where Db is the Jacobian matrix of the first measured 
variables yi=bi (x) . 

3 0 The terms null space, Jacobian matrix and inverse or 

pseudoinverse matrix, as applicable, have definitions which are 
familiar from the theory of matrix computations (see for 
example Gene H. Golub, Charles F. van Load: "Matrix 
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^ Computations", 3^^ Edition, Baltimore, London; The Johns 
Hopkins University Press; 1996) . 

In a further preferred form of embodiment of the invention, the 
5 first sensitivity variables calculated in the technical system 
are in each case the ratio of the change in accuracy of a 
selected parameter to the change in accuracy of a first 
measured variable, where the selected parameter is a selected 
state variable, which can "be determined via the first measured 
10 variables. The method is distinguished in this case by the fact 
that: 



- at least one of the selected parameters is a selected state 
variable, which can be determined via the first measured 

15 variables; 

- one or more of the first sensitivity variables O yjxi 
represents in each case the ratio of the change in accuracy 
Aoii^/xi = AE((xi - E(xi))^)/xi of the selected state variable 
xi to the change in accuracy Aojj^/yj = AE((yj - E(yj))^)/yj of 

20 a first measured variable yj; 

- the first sensitivity variables are determined from the 
following formula: 

$ = r^i,- 

25 YjXi o^u 

where rjj is the element in the 1^^ line and the j*^^ column of 
the matrix N-M'W. 



In a further form of embodiment, each of the first sensitivity 
30 variables represents the ratio of the change in the accuracy of 
a selected diagnostic variable to the change in the accuracy of 
a first measured variable, where the selected diagnostic 
variable can be determined via the first measured variables. In 
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this case, the method is distinguished by the fact that: 

- at least one of the selected parameters is a selected 
diagnostic variable, which can be determined via the first 

5 measured variables; 

- a matrix Dd is determined, this being the Jacobian matrix of 
the diagnostic variables di=di (x) ; 

- one or more of the first sensitivity variables ^ yj an 
represents in each case the ratio of the change in accuracy 

10 Aonn^/dn = AE ( ( dn - E(dn))^)/dn of the Selected diagnostic 
variable dn to the change in accuracy Aojj^/yj = AE((yj - 
E(yj))^)/yj of a first measured variable yj ; 

- the first sensitivity variables are determined by the 
'following formula: 

15 a'jj 

<^ = s^nj 

where s^j is the element in the n*^^ line and the j^^ column of 

20 

In a further preferred form of embodiment, one or more of the 
second sensitivity variables each represents the variance of a 
selected state variable when a second measured variable is . 
being added in, where the selected state variable can be 
25 determined via the first measured variables. The method is 
distinguished in this case by the fact that: 

- at least one of the selected parameters is a selected state 
variable which can be determined via the first measured 

30 variables; 

- one or more of the second sensitivity variables represents, 
in each case, the variance ak->xi^ of the selected state 
variable xi when a second measured variable, the value of 
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which is a state variable Xk with the variance Ck, is being 
added to the measurement park; 

- the second sensitivity variables are determined by the 
following formula: 

5 - (m^^-mi)' 

<^ ;c->xi = mi 'Ttii , 

where jni is the i^^ coliamn of the matrix M^-N. 

10 In a further form of embodiment of the invention, one or more 

of the second sensitivity variables represents in each case the 
variance of a selected diagnostic variable when a second 
measured variable is being added in, where the selected 
diagnostic variable can be determined via the first measured 

15 variables. Here, the method is distinguished by the fact that: 

- at least one of the selected parameters is a selected 
diagnostic variable which can be determined via the first 
measured variables; 

20 - a matrix Dd, which is the Jacobian matrix of the diagnostic 
variables di=di(x), is determined; 

- one or more of the second sensitivity variables represents, 
in each case, the variance ak->dn^ of the selected diagnostic 
variable dn when a second measured variable, the value of 

25 which is a state variable Xk and which has a variance Ok is 
being added to the measurement park; 

- the second sensitivity variables are determined by the 
following formula: 

30 (mj-qj' 

<^ k'>dn = Qn^'Qn " 
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where itii is the i^^ column of the matrix M^-lf, and qn is 
the n^^ column of the matrix and M^'N^'DdF. 

The case can now arise in which the selected parameter of the 
5 technical system is a state variable which cannot be determined 
via the first measured variables. In this case, the first step 
is to determine a second measured .variable, the value of which 
is a state variable, and which is to be added into the 
measurement park to enable the selected parameter to be 
10 uniquely determined. The method by which this case is taken 
into account is distinguished by the fact that: 

- at least one of the selected parameters is a selected state 
variable which cannot be determined via the first measured 

15 variables; 

- a matrix P, which is the orthogonal projection onto the null 
space of A, is determined; 

- a second measured variable is determined, the value of which 
is a state variable Xk, and which is to be added into the 

20 measurement park s.o that the selected state variable can be 
uniquely determined; 

- one of the second sensitivity variables represents the 
variance ak->xi^ of the selected state variable when the second 
measured variable x^ which has been determined, and which has 

25 the variance ojc, is being added to the measurement park; 

- the second sensitivity variable is determined by the 
following formula: 

WpV II IIpII II' 

30 o^k->xi = k + \\m i^k 11/ 

\\p4' II IIp/cII II 

with p = Pni, where ni is the 1^^ column of the matrix N^, and 
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mi is the i^^ column of the matrix M^-N^ and pk is the k^^ 
column of the matrix P-N^. 

r 

The further case can arise in which the selected parameter is a 
diagnostic variable which cannot be determined via the first 
measured variable. .In this case., the first step is to determine 
a second measured variable, the value of which is a state 
variable, and which is to be added into the measurement park. 
The method by which this case is taken into account is 
distinguished by the fact that: 

- at least one of the selected parameters is a selected 
diagnostic variable which cannot be determined via the first 
■measured variables; 

- a matrix Dd, which is the Jacobian matrix of the diagnostic 
variables di=di(x) , is determined; 

- a matrix P, which is the orthogonal projection onto the null^ 
space of A, is determined; 

- a second measured variable is determined, the value of which 
is a state variable Xk, and which is to be added into the 
measurement park so that the selected diagnostic variable can 
be uniquely determined; 

- one of the second sensitivity variables represents the 
variance aic->dn^ of the selected diagnostic variable dn when 
the second measured variable Xk which has been determined, 
and which has the variance a^c/ is being added into the 
measurement park; 

- the second sensitivity variable is determined by the 
f ol lowing formula : 



k->6n 




II 



Pk\ 



II 



2 



II 



II 



II 
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with p = Pcn, where is the n^^ coliimn of the matrix N^'Dcf, 
nik is the k^^ column of the matrix M^'N^ and p]^ is the k^^ 
column of the matrix P'N^. 

In the last two mentioned forms of embodiment, the 
determination of the second measured variable is preferably 
made by searching the matrix P'N*^ for the column for which p is 
linearly dependent on the vector which the column represents, 
and the index of this column then gives the second measured 
variable which is to be added to the measurement park. 

In a particularly preferred embodiment of the invention, in 
those forms of embodiment for which the variance of selected 
15 parameters is determined as the second sensitivity variable 
when second measured variables are being added in, 1% of the 
value of the second measured variable is taken as the standard 
deviation for the second measured variables which are to be 
added in. 

20 

It should be remarked at this point that the correctness of all 
the formulae used above can be demonstrated mathematically. 

In addition to the method described above, the invention also 
25 relates to a device for carrying out the method in accordance 
with the invention. Furthermore, the invention covers a 
computer program product which has a storage medium on which is 
stored a computer program, which can be executed on a computer 
and by which the method in accordance with the invention can be 
30 performed. 

Exemplary embodiments of the invention are explained and 
illustrated below by reference to the drawings 
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These show: 

Fig. 1 the schematic structure of a technical system which 
5 is analyzed by means of the method in accordance with 

the invention. 

Fig. 2 a processor unit for carrying out the method in 
accordance with the invention. 

The technical system shown in Fig. 1 relates to a heating 
system in a power station, with two heating surfaces 1 and 2 
connected one after the other, with a gas stream G and a water 
stream W flowing past these heating surfaces in opposite 
directions from each other. 

The technical system is characterized by the following state 
variables : 

20 mw,in mass flow of the water when it enters into the heating 
system; 

mw,out mass flow of the water when it leaves the heating 
system; 

25 

hw,in specif ic ' enthalpy of the water when it enters into the 
heating system; 

specific enthalpy of the water between the two heating 
surfaces 1 and 2 ; 



10 



15 



hw , middle 

30 



hw,out 



specific enthalpy of the water when it leaves the 
heating system; 
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ITIg 



mass flow of the gas in the heating system; 



he, in specific enthalpy of the gas when it enters into the 
5 heating system; 

ho, middle Specific enthalpy of the gas between the two heating 
surfaces 1 and 2 ; 

10 he, out specific enthalpy of the gas when it leaves the heating 
system. 

The state variables are the variables in a system of equations 
H(x)=Of which includes the following physical balance 
15 equations: 

Mass balance for the water in the heating system: 
20 Enthalpy balance at the first heating surface: 



Enthalpy balance at the second heating surface: 

fflff* r^Cjnlddle " ^,<mt) " in* (i^W, middle ~ t^W^in) =0. 

The following set operating points of the technical system are 
considered, where the values shown below for the state 
variables represent a solution to the above system of 
equations: 



25 









i^W, middle 








middle 




100 


100 


200 


300 


400 


50 


1000 


800 


600 



30 



Apart from the state variables identified above, the technical 
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system is further characterized by a diagnostic value. which, in 
the present case, represents the relative heat transfer of the 
gas which is flowing through. This heat transfer W can be 
described by the following formula: 

W = 

The following first measured variables are measured in the 
technical system with a standard deviation in each case of 1% 
relative to the setpoint value Concerned: 

Enthalpy flow for the water on entry into the heating system: 
Mass flow of the water on entry into the heating system: 
Enthalpy flow for the water on leaving the heating system: 
mass flow of the gas : 

Enthalpy flow of the gas on entry into the heating system: 

Using the formula quoted above for the relative heat transfer W 
with the setpoint values gives W = 0.4. 

Because the standard deviations used are 1%, the measurement of 
the relative heat transfer leads to a measured value of 0.4 
with a standard deviation of 0.0098. 
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In a first variant of the method in accordance with the 
invention, the sensitivity variable used for a first measured 
quantity is in each case the ratio of the change in accuracy of 
5 the diagnostic value W to the change in accuracy of the first 
measured value. 

This gives the following values: 

10 Sensitivity variable for the enthalpy flow of the water on 
entry into the heating system: 0.167. 

. Sensitivity variable for the mass flow of the water on entry 
into the heating system: 0.0. 
15 . 

Sensitivity variable for the. enthalpy flow of the water on 
leaving the heating system: 0.667. 

Sensitivity variable for the mass flow of the gas: 0.0 

20 

Sensitivity variable for the enthalpy flow of the gas on entry 
into the heating system: 0.167. 

It will be seen that the sensitivity variable for the enthalpy 
25 flow of the water on leaving the heating system has the 

greatest value. This means that a change in the accuracy of the 
measurement of the enthalpy .flow of the water on leaving the 
heating system has the greatest influence on the accuracy of 
the measurement of the relative heat transfer. It follows that 
30 an improvement in the measurement accuracy of the enthalpy flow 
of the water on leaving the heating system will be most 
effective in producing an improvement in the accuracy of the 
diagnostic value. By contrast, the measurements of the mass 
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flows have a sensitivity value of 0, and thus have no affects 
on the accuracy of the diagnostic value W. 

In a further form of embodiment of the method in accorciance 
with the invention, the variances of the relative heat transfer 
are calculated as the sensitivity variables, on the assiomption 
that in the technical system a state variable is being added in 
to the first measured variables as a second measured variable 
with a standard deviation of 1% relative to the setpoint value. 
Below are given the standard deviations (square roots of the 
variances) when individual state variables are being added in: 

Add in mvj^in' 0.0098; 

15 Add in mw^out- 0.0098; 

Add in hy^^in' 0.0095; 

Add in hf^^out: 0.0086; 

20 

Add in iogi 0.0098; 

Add in hcin- 0.0095; 

25 Add in he, out- 0.0062; 

The additi on of hw jxiiddie and h^ middle is not considered because 
these state variables are not uniquely determined by the 
measured variables. It can be seen that the introduction of the 
30 measurement of hc^out gives the smallest standard deviation for 
the relative heat transfer. As a consequence, the measurement 
of he, out is added into the measurement park for the first 
measurements . 



5 



10 
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In a further form of embodiment of the invention, consideration 
is now given to state variables which are not uniquely 
determined by the first measured variables of the technical 
5 system. In this case, these are the state variables h^, middle and 
hcmiddie- A first Step is now used to determine the measured 
variables which must be added in for the state variables 
hw.middie ^nd ha^middie to be uniquely determined. For this purpose, 
a calculation is performed in accordance with claim 11. 

10 

It turns out that a measurement of hw.middie or he, middle is in each 
case sufficient to determine the two state variables hw, middle 
and he, middle- Making the assumption of a standard deviation of 
1%' for the measurement of h^, middle or ho^middief as applicable, 
15 gives: 

- in the case of a measurement of hw.middiet a standard deviation 
for 

hw^middie of 3.0, and a standard deviation for ha, middle of 
17.34; 

20 

- in the case of a measurement of hG,middie» a standard deviation 
for hy^, middle of 9.06, and a standard deviation for ha, middle of 
8.0. 

25 From this it can be seen that for an exact determination of 

h]^, middle it is preferable to actually measure h^, middle/ whereas 
for an exact determination of ha, middle it is preferable to add 
also ho, middle ^s a measurement in the measurement park. 

30 The method described above permits a systematic and rapid 

search for measurement points by which the accuracy of selected 
state variables or diagnostic variables, as applicable, can be 
improved. It is thus no longer necessary to fall back on the 
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experience of engineers in order to decide which measured 
variables should be added into a measurement park, or which 
measurement accuracies should preferably be improved, as 
applicable. 

5 

Fig. 2 shows a processor unit PRZE for performing the method in 
accordance with the invention. The processor unit PRZE 
incorporates a processor CPU, a memory MEM, and an input / 
output interface lOS, which is used in various ways via an 

10 interface IPC: an output is shown visually on a monitor MON via 
a graphic interface, and/ or is output on a printer. Inputs are 
made via a mouse MAS or a keyboard TAST. The processor unit 
PRZE also provides a data bus BUS, which establishes the link 
between a memory MEM, the processor CPU and the input/output 

15 interface lOS. Further, additional components can be connected 
to the data bus, for example additional memory, data storage 
(hard disk) or scanners. 



20 



